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基于 mtDNA-COI 基因 序列 的 雷 氏 按 蚊 
分 子 群体 遗传 结构 研究 


Wo EJ. HEE 


(第 二 军医 大 学 病原 生物 学 教研 室 ， 上 海 200433) 


摘要 :【 目 的] 探讨 我 国 雷 氏 按 蚊 Anopheles lesteri 的 群体 差异 和 分 化 程度 。【 方法】 采用 PCR 方法 , 从 分 子 水 平 鉴别 了 
采 自 我 国 和 韩国 的 雷 氏 按 蚊 共 9 个 群体 139 个 样本 , 扩 增 其 线粒体 DNA. 细胞 色素 氧化 酶 亚 基 工 基因 , 并 进行 序列 测 
定 和 分 析 。[【 结果】 本 研究 共 获 得 49 个 单 倍 型 , 各 单 倍 型 呈 高 水 平 的 平行 演化 , 来 自 云 南 群 体 的 单 倍 型 显示 是 扩张 
的 源头 。 分 子 变异 等 级 分 析 (AMOVA) 的 计算 绪 果 显示 ,群体 内 变异 占 总 变异 的 比例 (64. 95% ) 大 于 群体 间 
(35.0596 ) ,Fsr 值 为 0.3504, 各 群体 间 出 现 遗 传 分 化 。Mantel 检验 结果 显示 基因 流水 平 与 地 理 距 离 呈 人 负 相 关 关 系 
(R^ 20.1322) , 群体 遗传 结构 符合 距离 隔离 模型 。【 结论] 雷 氏 按 蚊 韩国 和 辽宁 群体 与 其 他 分 布地 群体 间 差 异 大 , 已 
出 现 明显 分 化 , 我 国 其 他 分 布地 群体 间 的 遗传 差异 小 。 
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Molecular population genetic structure of Anopheles lesteri ( Diptera: 


Culicidae) based on mtDNA-COI gene sequences 

YANG Man-Ni, MA Ya-Jun' ( Department of Pathogen Biology, Second Military Medical University, 
Shanghai 200433 , China) 

Abstract. [ Aim] To study the population genetic variation and the diversity level of Anopheles lesteri from 
China. [Methods] The 139 specimen of 9 populations from China and South Korea were identified by PCR 
assay. Mitochondrial DNA cytochrome oxidase sununit I gene was amplified, and the mtDNA-COI gene 
sequences were sequenced and analyzed. [Results] In this study, 49 haplotypes were obtained, which were 
at a high level of homoplasy, and the haplotypes of Yunnan population showed to be the expansion origin. 
The results of the analysis of molecular variance ( AMOVA) showed that the percentage of variation within a 
population (64.95% ) was greater than that among populations (35.0596 ). The F, value was 0. 3504, and 
the genetic divergence among populations was significant. Mantel test results indicated that the level of gene 
flow was negatively correlated with geographic distances ( R^ 20. 1322). The population genetic structure fits 
the isolation-by-distance model. [Conclusion] The greatest difference is found between the populations of 
South Korea and Liaoning and other populations from China, and they have exhibited obvious differentiation. 
The level of variation among populations except those of Liaoning and South Korea is low. 
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ti Efi Anopheles lesteri Baisas and Hu, 1936 
隶属 赫 坎 按 蚊 种 团 (An，hyrcanus group) ， 如 同 按 蚊 
中 存在 的 其 他 隐 种 一 样 , 对 雷 氏 按 蚊 进 行 形态 分 类 
鉴别 非常 困难 , 因此 其 分 类 地 位 历经 多 次 修订 ( 软 
X&pt, 2008) 。 业 已 证 明 , 核糖 体 DNA 第 二 内 转录 
间隔 区 (ribosomal DNA second internal transcribe 
spacer, rDNA-ITS2 ) 序列 具有 客观 、 稳 定 的 种 间 差 


异 , 系 雷 氏 按 蚊 可 靠 \、 可 行 的 鉴别 特征 ( 马 雅 军 等 ， 
2005 ) 。 雷 氏 按 蚊 是 我 国 分 布 广 且 第 见 的 种 类 ,是 
tg DIE XR Z£ IR JN B5 38 92 Pr 8 AT C B E s 
S, 1997) , 但 在 不 同 地 理 分 布 区 其 生态 习性 和 传 
病 能 力 均 存在 一 定 差 异 ( 汤 林 华 , 2008) 。 有 鉴于 
此 , 其 不 同 地 理 群体 的 遗传 分 化 已 引起 有 关 学 者 的 
EM HES, 2002)。 
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分 子 群体 遗传 结构 是 近年 进展 较 快 的 热门 研究 
领域 ， 基 因 序 列 是 常用 的 分 子 标 志 ， 其 中 线粒体 
DNA 细胞 色 款 氧化 酶 亚 基 工 (mitochondrial DNA 
cytochrome oxidase subunit I, mtDNA-COI) 基因 序列 
被 用 于 多 种 按 蚊 的 群体 遗传 结构 分 析 , 成 为 具 重 要 
价值 的 分 子 标 志 之 一 (Walton et al., 2000; Chen et 
al., 2004; 王 冬 等 , 2007) 。 本 人 研究 在 雷 氏 按 蚊 分 子 
鉴别 的 基础 上 , 扩 增 和 分 析 mtDNA-COI 基因 序列 ， 
对 我 国 不 同 地 区 和 韩国 的 雷 氏 按 娄 群体 进行 遗传 多 
态 性 分 析 , 探讨 不 同 群 体 的 遗传 差异 和 分 化 现象 及 
其 影响 因素 , 为 败 媒 控制 提供 基础 资料 。 


1 材料 与 方法 


1. 1 按 蚊 来 源 

本 研究 样本 包括 2 个 实验 室 品 系 和 7 个 野生 群 
体 , 野生 样本 为 灯 诱 或 人 工 采 集 的 成 蚊 , 除 1 个 野 
生 和 群体 来 目 韩 国外 , 其 他 均 来 目 中国 ,， 群体 依据 采 
集 地 的 不 同行 政 省 份 命名 ( 表 1) 。 野 外 采集 的 成 蚊 
样本 , 依据 形态 特征 初步 鉴别 为 赫 坎 按 蚁 种 团 成 员 
PCR SEBESE, 1997) , 单 管 干燥 保存 市 回 实验 室 。 


表 1 本 研究 中 雷 氏 按 蚊 群 体 采集 信息 
Table 1 The collecting data of mosquito Anopheles lesteri populations in this study 


群体 采集 地 采集 时 间 ( 年 -月 ) 纬度 经 度 
Population Collecting site Collecting date ( Year-month ) Latitude ( N) Longitude ( E) 
海南 省 文昌 县 
HN* m 19°58’ 110°14 
Wenchang, Hainan Island 
GD 三 东 省 珠海 市 2007-10 22?17' 113?30' 
Zhuhai, Guangdong 
云南 省 盐 津 县 
YN Beim 2006-07 28?6' 104?13' 
Yanjing, Yunnan 
Vu JI ers TL TH NIS XE 
SC VERE TRAE TI 2006-07 28?1' - 30?14" 103?29' —104?34" 
Pujiang & Junlian, Sichuan 
湖北 省 广 水 市 和 随州 
HB 北 省 水 市 和 随州 市 2007-07 31?41' -31?52' 113?15' - 113?47' 
Guangshui & Suizhou, Hubei 
河南 省 
HEN TRIAL 2007 -07 32?29' 113?23' 
Tongbai, Henan 
TAE 
Js* NECS 32°54" 118*34' 
Xuyi, Jiangsu 
辽宁 省 东港 市 
LN í 2004-07 39°58” 123°50' 
Donggang, Liaoning 
Gang Hwe £ 
K SES Cang Hwe 8i 2005-07 37^6' 12790 
Gang Hwe Island, South Korea 
“实验 室 品系 ,其 数据 为 初始 群体 的 野外 采集 地 和 采集 时 间 ; 室 内 饲养 条 件 如 下 :温度 26 + 1 ,湿度 65% +5% 和 光 周 期 12L: 12D, " The 


laboratory colonies. The collecting site and date of their source populations from fields are listed here. They were bred at 26 +1°C and 65% +5% RH, 


under a 12L: 12D photoperiod in the laboratory. 


1.2. 雷 氏 按 蚊 的 分 子 鉴 别 

单 蚊 提取 基因 组 DNA 按 Collins 等 (1987 ) 的 方 
法 进行 , 将 所 提 的 DNA X F 100 pL 0X 2€ KAK 
中 , EEAICUKAB DUAE SES MYS rDNA-ITS2 序列 的 
差异 进行 PCR 分 于 鉴别 (Ma et al., 1998) ， 鉴 别 为 
雷 氏 按 蚁 的 样本 用 于 后 续 人 研究。 


1.3 mtDNA-COI 基因 扩 增 和 序列 测定 

以 雷 氏 按 蚊 基 因 组 DNA 为 模板 , 扩 增 mtDNA- 
COI 基因 片段 。 采 用 的 通用 引物 序列 依据 文献 
1994 ) 合成 : 正 向 5- 
GGTCAACAAATCATAAAGATATTGGC-3', Js |] 5”- 
TAAACTTCAGGGTGACCAAAAAATCA-3', f Ji pk 


( Folmer et al., 
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系 中 含 10 mmol/L Tris (pH 8.3) , 50 mmol/L KCI, 
1.9 mmol/L MgCL, 0. 1 mg/mL BSA, dNTPs 各 
0.1 mmol/L, DNA Taq 聚合 酶 2 U, 5%% 
0.1 pmol/L, 以 及 模板 1.5 Lo PCR 反应 条 件 为 : 
94C 2 min;94?C 30 s,45?C 30 s,72*C 30 s, 共 30 个 
1J83^;727C 8 min。 扩 增产 物 经 196 3 IR pase e Coo 
0.5 pg/mL 1R 45, 7, SE) 电泳 后 , $E I, 用 四 色 
荧光 标记 双 脱 氧 链 终止 法 测序 (ABI 3730 测序 仪 ， 
上 海 易 安生 物 科 技 有 限 公 司 ) 。 

1.4 统计 分 析 

测序 结果 以 Chromas 2. 13 进行 核对 确认 , 采用 
Clustal X 1. 8 进行 序列 比 对 ,， MEGA 2. 1 软件 包 统 
计 序 列 特征 , 单 倍 型 之 间 的 95% 简约 家 系 网 络 关系 
由 TCS 1.21 构建 。 

将 比 对 的 序列 数据 转换 为 ARLEQUIN 3.0 可 识 
别 的 文件 后 , 统计 和 计算 各 群体 的 序列 碱 基 分 化 参 
数 , 包括 置换 位 点 数 (5)、 分 离 位 点 数 (9s) 和 错 配 平 
均 数 (09. ) (Schneider et al., 2000) ,以 及 成 对 群体 间 
的 遗传 距离 。 依 据 遗 传 距离 由 MEGA 2.1 构建 群体 
间 的 聚 类 关系 。 群 体 遗 传 结构 的 参数 由 ARLEQUIN 
的 分 子 变 异 等 级 分 析 (analysis of molecular variance, 
AMOVA ) 模 块 计算 , 包括 :群体 间 的 分 化 参数 Pr( F- 
statistics) 基因 流水 平 Nm(Nm 21 - F4/4F,) 群体 
内 和 群体 间 的 变异 成 分 比例 。 由 ARLEQUIN 相应 的 
模块 进行 错 配 分 析 (mismatch analysis) ， 中 性 检验 
( neutrality test) H1 Tajima’ s } D FU Fu’ s 的 Fs 衡量 。 
Genepop 1.2 进行 Mantel 检验 , 分 析 地 理 距 离 与 群体 
分 化 程度 的 相关 性 。 

祖先 群体 的 规模 (0) 由 以 下 公式 计算 :0 =2Nu， 
其 中 N 代表 群体 内 有 效 肉 蚊 数 , u 代表 每 代 的 突变 
K, u =/ 由 (由 是 每 位 点 每 代 的 突变 率 , 左 是 分 析 序 列 
的 碱 基数 ) 。 另 外 , 应 用 7 -2ut 可 估算 群体 扩张 的 近 
似 世 代数 ( (Rogers and Harpending, 1992) , t 是 根 
据 广义 非 线性 最 小 二 乘法 估计 出 的 群体 扩张 参数 。 


2 结果 


2.1 雷 氏 按 蚊 的 mtDNA-CO I 基因 序列 特征 

本 研究 扩 增 获得 139 条 雷 氏 按 蚊 的 mtDNA-COIT 
基因 序列 , 对 比 后 采用 其 中 长 度 为 660 bp 的 片段 进 
行 后 续 分 析 , 该 片段 的 GC 含量 平均 为 31.4% , fà 
基 置 换 位 点 52 个 ,其 中 34 处 为 转换 (65. 38% ) ， 
15 处 为 颠 换 (28. 85% ) , 3 个 超 变 位 点 位 于 变异 位 
点 的 65、124 和 556 处 , 无 插入 和 缺失 ( 表 2), 共 49 
个 单 倍 型 。 


群体 中 所 有 样本 的 单 倍 型 分 化 水 平 与 核 苷 酸 序 
列 的 差异 程度 密切 相关 , 核 苷 酸 的 分 化 程度 依据 8。 
和 0, 衡量 , 野生 群体 的 0,(2. 18 +1.49 - 6.49 + 
2.89) $10,(2.00 +1.68 ~5.01 £2. 77) H] 3E 
验 室 品系 ( 表 3)。 
2.2. 雷 氏 按 数 mtDNA-COI 单 倍 型 间 的 关系 及 其 
在 群体 中 的 分 布 

本 研究 结果 显示 雷 氏 按 蚊 的 mtDNA -COI. 基因 
序列 存在 49 个 单 倍 型 ， 所 后 比例 为 35. 25% 
(49/139) 。 其 中 在 群体 间 共 享 的 单 倍 型 为 7 个 , 分 
别 为 :hp2, h3, h4, h6, h8, h9 和 h45, 占 整 个 单 倍 
型 的 14.29%( 表 3) 。 单 倍 类 型 在 群体 中 分 布 的 个 
体 超过 10 个 的 有 h6 (20. 1496, 28/139) ) 和 h9 
(10.79% , 15/139) , 仅 1 个 个 体 的 占 63.27% (31/ 
139)。 雷 氏 按 蚊 单 倍 型 家 系 网 络 图 显示 (图 1), 各 
单 倍 型 呈 高 水 平 的 平行 演化 , 提示 雷 氏 按 蚊 群体 在 
历史 上 曾 发 生 多 次 局 部 扩张 。 来 自 云 南 群 体 的 单 倍 
型 (h2, h4, h6) 显 示 是 扩张 的 源头 , 经 过 1 ~2 PR 
变形 成 广东 和 江苏 群体 中 的 单 倍 型 ;经 过 3 步 以 上 
的 突变 形成 辽宁 和 韩国 群体 的 单 倍 型 , 在 此 之 前 的 
单 倍 型 混杂 分 布 于 河南 、 四 川 和 湖北 群体 。 图 中 显 
示 未 检测 到 的 单 倍 型 44 个 , h2 作为 外 群 的 最 大 可 
能 值 是 0.082 。 





图 1 雷 氏 按 蚁 群体 的 mtDNA-COI 基因 单 倍 型 家 系 网 络 图 
Fig. 1 Genealogical relationship among mtDNA-COI 
gene haplotypes of Anopheles lesteri populations 
圆圈 代表 单 倍 型 , 圈 内 的 数字 代表 单 倍 型 的 类 型 , 之 间 的 连 线 代 表 
一 个 突变 步 又 , 空白 圆圈 是 未 检测 到 的 单 倍 型 。The circle indicates 
one of haplotype and the number in circle represents the type of 
haplotype. Each mutation step is shown with a short line. The empty 


circles indicate haplotypes that were not observed. 
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R2 雷 氏 按 旷 群体 各 mtDNA-COI 单 倍 型 碱 基 差 异 
Table 2 Variable bases of the mtDNA-COI gene for the haplotypes of Anopheles lesteri populations 
变异 位 后 ， 

单 倍 型 Variable position EU 
Haplotypes 11112222223333333333444444555555555556666 666 GenBank 

3366677700240012661334567990123569005557777800001122 accession no. 

2405603949453816249470544173240799571363456714670156 
COI hl cACITrITTCTTT[AIT TATCTCGCTCTGGATCAGAAGTGTGTÍ[A[TAACTTTACAGCC EU699001 
COI m ...[]..... a ellen LL... C....... EU699002 
COT h3 E- n LN G...... [ C..... T. EU699003 
COI h4 BO... EJ oaa G...... IL... C....... EU699004 
CO h5 E- [LT O lllo LLE.-.-e- C....... EU699005 
COI h6 S a a n Ti G...... H..... C. . T. EU699006 
COI h7 ELb....- FL.-.......... Asse. IL... C. . T. EU699007 
COI h8 -El..c..RLHbe........- Teens A- .A...[l}..... C. . T- EU699008 
COI h9 ...[l1Cc SS a A..... A. .A...H.... C....... EU699009 
COI hl0 ...[] ..... LL...... A...AA.............. a C....... EU6990010 
COT hll .器 -T.-H LE... CC...... EU6990011 
COI hl2 ...[]..... LIL...... Asse A. cH .... C..... T. EU6990012 
COI hl4 A. .[1..... [KG [] .... C....... EU6990014 
col hi5 ...[]....- Fi... A........ [l -.... C--T-- T- EU6990015 
COI hi6... Eb. T...El ooa ILH....- C....... EU6990016 
COI hl7 ...[-..-... a G...... H-T. AC- ATG EU6990017 
COI m0 ...[1]1..... LEC.G......... llle [jj.-.... C....... EU6990020 
COI mi ...[]l..... FL...... A. ees G-...... El-.... C....... EU6990021 
COI h22 ...Ll..... [IT GA. .... El-.-...-. C....... EU6990022 
col h23 A..[]..... E GA... .. LL... C....... EU6990023 
COI h28 . . . [A]A. "E ci G-. ..... El..-....- C..T. EU6990028 
COI h29 . .TC --- A- EE Too nnnm G... s. EH- C. . T. EU6990029 
COI h30 ...|cl ..... a G...... Ej..... C.. T. EU6990030 
COI n3 ...Ll..... LIE... Teeeeeeeee T- G- El....- C..T.... EU6990031 
col np32...L[ 1] ..... S s G...... Ll.... C.. T.. T. EU6990032 
COI n3 ...[E]....- Ll-.G.-.-..-..-.-..... C. .G- C... ÀBl...- C....... EU6990033 
COI h34 ...LEl..... E e A-.-..-.-... Lee C....... EU6990034 
col b35... [}p..... LIl.G.........2.....- A........ H- C....... EU6990035 
COI n36 --- [LT IT LlcG...c.. T.. T- EU6990036 
COI h37 --. El -..-...- LF. EH- C-.T-- T- EU6990037 
COI n38 ...Ll]1..... EF.-...-... aaa aaa A........ El.... C....... EU6990038 
COI n39 ... [E]... Eb... T eene G-...... HQ- C-....... EU6990039 
COI n40.G.Ll..... El...T.............. TG...... [E]... C....... EU6990040 
COI p41 - G-E}... FL...TC...-.......... TG-..-.... H- C....... EU6990041 
COI h42 .-- [lc G- eee Hb- C-- T- EU6990042 
COI h43 ...L[1..... orr G-...... Ll.... C-....... EU6990043 
COI h44 --.[l..... ILI... Teeeeeee- A...... G - c- H --A-C- G... EU6990044 
COI h45... E]... Eb... T....... A........ Ll... C.. T. EU6990045 
COI há6 .- E]... Ll...... A SA A- G...... Ll.... C..T.... EU6990046 
COI h47 .-.-.[|...... Hi-...-.- AT- -- A- GA- G- - A- EL]... -.- C.. T.. T. EU6990047 
COI h48 ... LE] ..... [L .---T--A----A---T-A--A----- Fl.... C..T.. EU6990048 
COI h49 ...[]....- LHLb..........-. A..... A... A. SS Ll. C.. TG. EU6990049 
COI hs0 --. El c... E] eese A.......... A.G...... BG. . ... C.. T. EU6990050 
col ht .CC. Ep... T....... A.... A... Bl... C.. T. EU6990051 
COI hs -- L1 C... EE eese As Asl A... AE eese Ces EU6990052 
COI hs --.-.Elc.... Bb... A.-......... A. G. G..Bl...- C.. T- EU6990053 
COI hs4 --.. El c.... B... A.......... A.G...... 回 -.... C..T.... EU6990054 
COI h55 .--.[}..... LL....-- AT...A..... A.G.. A... Ll]... C-. T.. T. EU6990055 
COI h56 -.-. El ..C.. El... sss OPE A. G...... LL...- C.. T. EU6990056 





数字 表示 序列 比 对 的 位 置 ,“. ”表示 与 单 倍 型 COL_hl 的 碱 基 相 同 , 方 框 为 超 变 位 点 。Numbers indicate the alignment positions, the small dots 
indicate identity with the COI_h1 haplotype, while the frames indicate the hypermutation positions. 


2.3 


雷 氏 按 蚊 的 群体 遗传 结构 
对 7 个 野生 群体 进行 遗传 结构 分 析 , 结果 显示 
成 对 群体 间 的 基因 流 范围 为 0. 136( HB-K) ~ 6. 893 
(YN-HEN) (4X 3) , 总 体 趋势 是 基因 流水 平 随 者 地 
理 距 离 的 增 大 而 减 小 。 对 群体 间 的 地 理 距 离 与 基因 


流 的 相关 性 进行 回归 的 Mantel 检验 , 结果 显示 基因 
流 与 地 理 距 离 呈 负 相 关 关 系 (R =0. 1322, P = 
0.02951) (El 2) , 群体 遗传 分 布 格局 符合 距离 隔离 


模型 (isolation-by-distance model) 。 


1004 


群体 


Population 


CD 


YN 


SC 


HB 


HEN 


JS 


LN 


K 


样本 数 
Sample 


size 


10 


24 


11 


24 


35 


4 


单 倍 型 (n) 
Haplotypes 


h9 * (4) 


h39(9) , h40(4) , 
h4l(1), h42(2), 
h43(1) , h44(1), 


h45 * (2) 


h1(3), h2* (2), 
I8 * (1), h4" (1), 
h5(1), h6* (1), 


h7(1) 


h33(15) , h34(1) , 
h35(3) , h36(1), 
h37 (1) , h38(3) 


h4* (2), h20 (2), 
h21(1), h22(1), 


h23(1) 


I2 * (2) ,h3 (1), 
h10(2) , h11(1), 
h12(1) , h14(1), 
hl15(1) h16(1), 


h17(1) 


h6* (15), h28(1), 
h29(1), h30(1), 
h31(1), h32(5) 


h2* (2), h6* (12), 
h8* (1), hb9* (3), 
h45 * (4), h46(1), 
h47(1) , h48(4) , 
h49(1) , h50(1), 
h51(1) , h52(1), 
h53(1) , h54(1), 


h55(1) 


h8* (3), h56(1) 
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Table 3 Summary of Anopheles lesteri populations assayed 


N 序列 间 平 均 错 配 群体 内 和 群体 间 的 遗传 距离 ” 
分 离 位 点 数 (b， 
置换 位 点 数 (S) 7 ) 平均 数 ( 0 ) Average Genetic distances within 


Number of segregating 


Number of D pairwise number and between populations 
， nucleotide sites per | 
segregating sites of differences between 


sequence (+ SD) 


sequences (* SD) GD YN SC HB HEN LN K 
0 0.00 x0. 00 0.00 +0. 00 - - - - - - 一 
14 3.95 +1.66 2.82 +1.73 2.000 0. 966 3.329 0.923 1.444 1.245 4.390 
7 2.47 +1.31 2.09 +1.44 0.268 2.088 1.920 0.346 0.137 1.183 3.606 
12 3.21 +1.36 2.94 +1.78 0. 586 0. 412 2.938 1.990 1.707 3.040 5.156 
6 2.45 +1.41 2.48 +1.74 0. 249 0. 137 0.409 2.476 0.336 1.577 4.262 
19 6.49 x 2. 89 4.62 x2. T] 0. 301 0.035 0. 337 0.068 4.618 1.276 3.782 
8 1. 88 +0.90 1.00 x0. 78 - - - - 一 一 一 
20 4.86 +1.77 5.01 +2.77 0. 224 0. 197 0. 420 0.236 0.204 5.008 1.976 
4 2.18 +1.49 2.00 +1.68 0.614 0. 635 0.641 0.647 0.470 0.254 2. 000 


“共享 的 单 倍 型 , 括号 内 的 数字 代表 该 单 倍 型 在 群体 中 的 个 体 数量 。” 上 三 角 为 群体 间 校 准 的 核酸 差异 平均 数 ; 对 角 线 为 群体 内 差异 的 平均 
数 ( 黑体 字 ) ;下 三 角 为 基于 群体 间 和 遗 传 距离 计算 的 Fst f," Haplotypes occur in more than one population and the number in parentheses indicates 


the frequency of the haplotype in the population. "' The average number of pairwise differences within a population and corrected average pairwise 


difference between populations are, respectively, on and above the diagonal; genetic distance between populations calculated as pairwise Fst values is 


below the diagonal. 
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图 2 基因 流 与 地 理 距离 间 的 相关 性 分 析 
Fig. 2 Correlations between gene flow level and 


the geographic distance among populations 


分 子 变异 等 级 分 析 的 计算 结果 显示 , 群体 内 和 
群体 间 均 有 较 高 的 遗传 变异 , 对 总 变异 的 贡献 前 者 
多 于 后 者 ， 和 群体 间 发 生 显著 遗传 分 化 (Pr 值 为 
0.3504, 表 4)。 依 据 遗 传 距离 构建 雷 氏 按 蚊 各 群 
体 的 聚 类 关系 (图 3) ,结果 显示 韩国 和 辽宁 群体 与 
我 国 其 他 分 布地 的 群体 位 于 不 同 的 分 支 , 其 间 差 异 
K, 出 现 显 著 遗 传 分 化 , 而 其 他 分 布地 群体 间 的 遗 
传 差异 小 。 


表 4 雷 氏 按 蚊 群体 遗传 变异 的 分 子 变 异 等 级 分 析 结 果 
Table 4 AMOVA analysis of genetic variation in Anopheles lesteri populations 


遗传 变异 来 源 自由 度 
Source of variation Degree of freedom 
群体 间 6 
Among populations 
群体 内 104 
Within a population 
总 计 Total 110 


固定 指数 Fixation indiex; Fsr =0. 3505. 





0.5 


图 3 雷 氏 按 蚊 各 群体 间 的 聚 类 关系 (NJ 法 ) 


Fig. 3 The clustering relationship among Anopheles lesteri 


populations based on pairwise genetic distance ( NJ method) 


应 用 Arlequin 软件 中 的 中 性 检验 模块 ， 计 算 
Tajima’ s BJ D FI Fu’ s 的 Fs (&, 结果 显示 , 除了 湖 
北 和 辽宁 群体 的 D 为 正 数 外 , 其 他 的 均 为 负数 ( 表 
5), 提示 群体 内 存在 大 量 的 低频 突变 , 群体 正在 处 
于 扩张 阶段 。 歧 点 分 布 图 显示 群体 单 倍 型 碱 基 差 异 
分 布 为 明显 的 单 峰 (图 4), 提示 为 群体 符合 扩张 模 
型 , 方差 检验 结果 显示 检测 结果 与 期 望 值 的 符合 程 
度 无 统计 学 意义 (SSD = 0. 0018, P = 0. 60) 
( Schneider and Excoffier, 1999) 。 依 据 公 式 计 算 7 = 
0.00 ~ 92. 20, 祖先 群体 规模 0, = 0. 00 ~ 0. 39， 
0 = 0 ~99 999, 参照 文献 采用 mtDNA-COI 基因 每 
位 点 每 代 的 突变 率 1 x 10 ”为 雷 氏 按 蚊 的 进化 速率 
( Powell et al., 1986) , 估算 群体 内 有 效 肉 蚊 数 约 为 
0.20 x10 , 若 以 按 蚊 在 1 个 月 内 发 生 2 个 世代 计 
算 , 则 全 局 性 的 群体 扩张 大 约 发 生 在 582 000 年 前 。 


遗传 变异 元 素 占 总 变异 百分比 (%) 


Variance components Percentage of variation 


0.9789 35.05 
1.8140 64.95 
2. 7929 100 


RS 雷 氏 按 蚊 群体 中 性 检验 数值 
Table 5 Values of neutrality test for Anopheles 


lesteri populations 


指数 Index GD YN SC HB HEN LN K 


Tajima's D -1.04 -0.66 -0.29 0.06 -1.30 0.11  -0.78 


Fu's Fs -25.89-10.31 -26.38 -5.02 -7.41 -25.53 -1.87 


--*-- 观 察 值 Observed value 
A —a—Jdi[ Simulated value 





出 现 频率 Frequency 
B 
e 
© 


9 10 11 12 13 14 15 
MEAM Nucleotide difference 


图 4 雷 氏 按 蚊 群 体 mtDNA-COI 单 倍 型 的 政 点 分 布 图 
Fig. 4 Mismatch distribution of the observed mtDNA-COI 


1 2 3 4 5 6 7 8 


haplotypes of Anopheles lesteri population 
3 讨论 
雷 氏 按 蚊 在 我 国 日 前 记录 的 范围 包括 100°E 以 
ZR, 19°N ~ 42°N 的 18 个 省 (市 /区 ) CUR AB, 


2008) 。 本 研究 对 雷 氏 按 蚊 群 体 的 收集 策略 是 尽量 
ARDRE, 采集 的 按 蚁 样本 经 分 子 鉴别 确认 
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为 雷 氏 按 蚊 后 进行 群体 遗传 结构 研究 。 各 采集 地 按 
蚊 的 分 子 鉴别 结果 显示 , 除 辽 宁 和 广东 珠海 横 容 岛 
外 ,其 他 地 点 的 雷 氏 按 蚊 群 体 密度 较 低 或 无 雷 氏 按 
蚊 ( 如 山东 ` 贵 州 \ 广 西江 苏 和 福建 等 ) ， 且 我 们 的 
研究 ( 男 文 发 表 ) 还 发 现 有 些 采 集 地 雷 氏 按 蚊 及 其 
隐 种 的 组 成 和 比例 与 以 往 依据 形态 鉴别 的 结果 不 同 
(FLSE, 1993; AIRE, 2008), 值得 进一步 探 
讨 。 虽 然 如 此 , 本 研究 还 是 获得 了 分 布 于 我 国 东 
北西 南华 南 和 长 江淮 河流 域 等 的 雷 氏 按 蚊 野生 和 群 
体 , 样本 应 具有 较 强 的 代表 性 。 

本 研究 应 用 mtDNA-COI 基因 序列 分 析 雷 氏 按 
蚊 群 体 的 遗传 变异 和 群体 结构 ,结果 显示 该 基因 序 
列 突变 位 点 分 布 均匀 ,野生 群体 的 核 苷 酸 分 化 程度 
大 于 实验 室 品系 , 进一步 证 明 mtDNA-COI 基因 是 
按 蚊 群体 遗传 分 析 的 理想 分 子 标志 (Walton et al., 
2000; Chen et al., 2004; 王 冬 等 ,2007) 。 

漳 祖 理论 认为 ,从 群体 中 所 抽样 本 沿 时 间 往 回 
追溯 分 析 , 可 以 推 采 这 些 样 本 的 最 近 共 同 祖先 (the 
most recent common ancestor ) 产生 以 后 的 时 间 段 内 
所 发 生 的 事件 , 即 从 现存 群体 中 抽样 可 推论 群体 的 
历史 。 本 研究 雷 氏 按 蚊 个 体 的 mtDNA-COI 序列 高 
度 相 似 的 现象 , 用 漳 祖 理论 可 很 好 地 解释 , 即 它们 
分 享 共 同 祖先 的 mtDNA-COI 序列 ;祖先 群体 分 化 
后 , 核 音 酸 序 列 在 漫长 的 时 间 内 不 断 发 生 突 变 , 并 
随 着 气候 和 地 理 因 子 的 变化 ,以 及 群体 扩张 和 迁移 
等 , 形成 目前 的 群体 遗传 结构 格局 。 综 合 本 研究 结 
果 提 示 雷 氏 按 蚊 可 能 曾经 历 过 近期 群体 扩张 事件 ， 
且 目 前 仍 处 于 群体 扩张 阶段 , 故 mtDNA-COI 基因 
序列 的 核 苷 酸 变 异 主 要 发 生 在 于 群体 内 个 体 间 。 

固定 指数 Fsr 通 常 是 衡量 不 同 群体 间 分 化 程度 的 
指标 , 基因 流 Nm 则 可 衡量 群体 间 基 因 交 流 的 水 平 ， 
本 研究 结果 提示 总 体 的 趋势 是 各 群体 之 间 的 基因 交 
流 受 到 一 定 程 度 的 限制 , 群体 间 地 理 距离 愈 远 , 基因 
流水 平 愈 低 。 理 论 上 , 基因 流水 平 小 于 1 则 难以 抵 
御 遗 传 漂 变 造成 的 群体 间 遗 传 分 化 (Slatkin，1987)， 
本 人 研究 除 云南 与 湖北 \ 河南 .辽宁 群体 之 间 ， 以 及 河 
南 与 湖北 群体 之 间 的 基因 流水 平 大 于 1 外 , 其 他 的 
群体 之 间 的 基因 流水 平均 较 低 ( 小 于 1), 各 群体 的 
聚 类 关系 可 看 出 咒 国 和 辽宁 群体 与 我 国 其 他 分 布地 
群体 间 差 异 较 大 ,已 出 现 明 显 分 化 ,这 与 应 用 
RAPD 分 析 的 结果 一 致 ( 马 雅 匣 等 , 2002 ) 。 

群体 遗传 结构 形成 的 影响 因素 应 与 分 布地 的 气 
修 和 地 理 特征 , 以 及 蚊虫 的 生态 习性 等 结合 后 综合 
分 析 。 众 所 周知 , 我 国 云 南 省 地 理 地 貌 复 杂 ， 是 从 


南 喜 玛 拉 雅 山 脉 到 东亚 、 从 东南 亚 热 寓 到 中 国 的 亚 
热带 气候 的 过 渡 区 ,同时 又 是 从 内 瓦 纳 古 大 陆 和 欧 
亚 板 块 分 离 后 的 印度 板块 和 缅甸 板块 的 结合 处 
( Audley-Charles, 1987) , 被 称 为 著名 的 生物 多 样 性 
的 中 心 ( 董 学 书 等 , 2004 ;Zhu et al., 2004) 。 在 冰川 
时 期 , 云南 可 能 保存 了 大 量 中 生 代 蚊 忠 的 栖 居 地 ， 
致使 体型 弱小 的 蚊虫 能 够 生存 下 来 ， 当 旱季 来 临 
时 , 蚊虫 随即 扩散 。 假 设 雷 氏 按 蚊 的 祖先 群体 来 目 
云南 ,其 他 群体 均 为 冰川 后 期 从 云南 扩散 而 至 , 结 
合 本 研究 的 单 倍 型 家 系 网 络 图 提示 , 雷 氏 按 蚊 从 云 
南 癌 北 、 东 迁移 和 扩散 ,在 途中 定 殖 并 出 现 群 体 局 
部 扩张 。 从 云南 回 东 北 迁 移 和 扩散 的 线路 中 经 历 的 
湖北 、 河 南 和 四 川 群体 之 间 虽 然 存 在 一 定 的 地 理 屏 
障 , 但 未 影响 其 间 的 基因 交流 ， 其间 遗传 差异 较 
小 ,尚未 形成 明显 的 遗传 分 化 ; 定 殖 在 辽宁 和 韩国 
的 群体 是 经 过 多 步 突 变形 成 , 与 祖先 群体 差异 最 
大 。 当 然 知 要 证 实 上 述 雷 氏 按 蚊 群体 的 迁移 和 扩散 
路 线 ， 尚 需 增加 云南 以 南 的 样本 进一步 证 实 。 雷 氏 
按 蚊 为 我 国 古 北 界 和 东洋 界 的 广 布 种 , 各 分 布 区 气 
候 条 件 不 同 , 且 人 区 生地 和 嗜 血 习性 也 存在 差异 Jl 
我 国 南部 和 中 部 的 雷 氏 按 蚊 主要 嗜 吸 人 血 , 而 在 北 
部 辽宁 则 偏 吸 动物 血 ( 汤 林 华 ,2008) ;加 之 由 于 大 
量 杀 虫 剂 的 使 用 , 导致 雷 氏 按 蚊 的 群体 数量 小 , dE 
因 交 流水 平 不 高 。 综 合 历史 和 现状 分 析 ， 上 述 因素 
可 能 是 影响 雷 氏 按 蚊 的 群体 遗传 结构 的 关键 。 

据 重 新 修订 的 记录 , 雷 氏 按 蚊 分 布 于 中 国 \ 日 
本 关岛 、 转 国 和 菲律宾 等 地 (Rueda et al., 2005) , 
从 分 子 特征 上 也 已 证 明 均 为 同 种 ( Wilkerson et al., 
2003; Hwang et al., 2006) , 但 就 医学 重要 性 而 言 ， 
雷 氏 按 蚊 在 关岛 和 非 律 宾 的 町 咬 习性 不 明 , 未 被 当 
作 传 症 媒 介 ， 而 在 我 国 中 部 却 具 有 较 高 的 媒介 效能 
( 陆 宝 甩 等 ,1997 ; 汤 林 华 , 2008) , 在 罪 国 和 日 本 也 
具有 一 定 的 传 症 作 用 (Tanaka et al., 1979; Shin et 
al., 2002) 。 因 此 , 加 强国 际 合作 , 应 用 更 多 的 分 子 
标志 ,系统 深入 地 研究 雷 氏 按 蚊 的 群体 遗传 , A RE 
找到 控制 传 症 作 用 相关 表 型 的 基因 , 将 为 媒介 控制 
措施 的 制定 起 到 重要 的 指导 作用 。 
致谢 承蒙 辽宁 省 东港 市 疾病 控制 中 心 陈 哲 、 江 苏 
省 寄生 来 病 防治 所 周 华 云 四川 省 疾病 控制 中 心 雷 
心田 和 韩国 Inha 大 学 Gi-Sik Min 教授 赠送 样本 , x 
南 省 寄生 来 病 防治 所 顾 云 安 、 湖 北 省 疾病 控制 中 心 
黄 光 全 和 陈 国 英 、 河 南 省 疾病 控制 中 心 夺 鞍 、 吐 州 省 
从 江上 县 疾病 控制 中 心 石 跃 华 、 广 东 省 珠海 市 疾病 控 
制 中 心 黄 利 群 及 第 二 军医 大 学 攀 勇 和 吴 静 等 协助 采 
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集 样本 ， 在 此 一 并 表示 诚挚 的 谢意 ! 
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